Stability of f{R) black holes 



fl, and Edwin J. Sonfl 



Yun Soo MyungxIjTaeyoon Moont], and Edwin J. Son 

Institute of Basic Sciences and School of Computer Aided Science, Inje University 

Gimhae 621-749, Korea 
^ Center for Quantum Space-time, Sogang University, Seoul, 121-742, Korea 

Abstract 

We investigate the stabihty of f{R) (Schwarzschild) black hole obtained from the 
f{R) gravity. It is difficult to carry out the perturbation analysis around the black 
hole because the linearized Einstein equation is fourth order in f{R) gravity. In order 
to resolve this difficulty, we transform f{R) gravity into the scalar-tensor theory by 
introducing two auxiliary scalars. In this case, the linearized curvature scalar becomes 
a scalaron, showing that all linearized equations are second order, which are the same 
equations for the massive Brans-Dicke theory. It turns out that the f{R) black hole 
is stable against the external perturbations if the scalaron does not have a tachyonic 
mass. 
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1 Introduction 



Modified gravity theories, f{R) gravities [H El [3] have much attentions as one of strong 
candidates for explaining the current accelerating universe f{R) gravities can be con- 
sidered as Einstein gravity (massless graviton) with an additional scalar. For example, it 
was shown that the metric-/(i?) gravity is equivalent to the wbd = Brans-Dicke (BD) 
theory with the potential [5]. Although the equivalence principle test in the solar system 
imposes a strong constraint on f{R) gravities, they may not be automatically ruled out if 
the Chameleon mechanism is introduced to resolve it. It was shown that the equivalence 
principle test allows f{R) gravity models that are indistinguishable from the ACDM model 
in the background universe evolution [6]. However, this does not necessarily imply that 
there is no difference in the dynamics of perturbations [7]. 

In order for f{R) gravities to be acceptable, they must obey certain minimal require- 
ments for theoretical viability [H [S] . Three important requirements are included: (i) they 
possess the correct cosmological dynamics, (ii) they are free from instabilities (tachyon) 
and ghosts [8l|9l[6], (iii) they attain the correct Newtonian and post-Newtonian limits. 

On the other hand, the Schwarzschild-de Sitter black hole was obtained for a positively 
constant curvature scalar in [7] and other black hole solution was recently found for a 
non-constant curvature scalar [10]. A black hole solution was obtained from f{R) gravities 
by requiring the negative constant curvature scalar R = R [H]. If 1 -|- f'{R) > 0, this 
black hole is similar to the Schwarzschild-AdS (SAdS) black hole. In order to obtain the 
constant curvature black hole solution from f{R) gravity coupled to the matter, the trace 
of its stress-energy tensor T^,y should be zero. Hence, two known matter fields are the 
Maxwell [I2] and Yang-Mills fields [I3]. 

All black hole solution must pass the stability test. A black hole solution should be 
stable against the external perturbations because it stands as the physically realistic object. 
There are two ways to achieve the stability of a black hole: one is the full stability by 
considering odd and even perturbations [H] and the other is the restricted stability by 
taking into account the spherically symmetric perturbations for simplicity [15]. The latter 
is not enough to guarantee the full stability and thus, it must be supported by the black 
hole thermodynamics (heat capacity) [T^. The basic idea is to decouple the linearized 
(perturbed) equations and then, manage to arrive at the second order Schrodinger-type 
equations for the physical field with the potential. If all potentials are positive for whole 
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range outside the event horizon, the black hole under the consideration is stable. Studies of 
stability of Kerr black hole are not as straightforward [T7] , because it is axially symmetric 
black hole and thus, the decoupling process seems to be complicated. However, this method 
is not suitable for f{R) black holes because f{R) gravity contains fourth order derivatives 
in the linearized equations [HI [19]. In this case, the requirement (ii) will play an important 
role in testing the stability of f{R) black holes. 

In this work, we investigate the stability of f{R) (Schwarzschild) black hole. We trans- 
form /(-R) gravity into the scalar-tensor theory to eliminate fourth order derivative terms by 
introducing two auxiliary scalars. Then, the linearized curvature scalar becomes a scalaron, 
indicating that all linearized equations are second order. Interestingly, they are exactly the 
same equations for the massive Brans-Dicke theory. Using the stability analysis of black 
hole in the massive Brans-Dicke theory, we show clearly that the f{R) black hole is stable 
against the external perturbations if the scalaron does not have a tachyonic mass. 

2 Perturbation of f(R) black holes 

Let us first consider f{R) gravity without any matter fields whose action is given by 



where ' denotes the differentiation with respect to its argument. It is well-known that 
Eq. (12.21) has a solution with constant curvature scalar R = R. In this case, Eq. (12. 2p can 
be written as 




(2.1) 



where k"^ 



SttG. The Einstein equation takes the form 




(2.2) 



RfiufiR) 



^g,J{R) = 0, 



(2.3) 



and thus, the trace of (12. 3p becomes 



RfiR) 



2f{R) = 0. 



(2.4) 



Note that the above equation determines the constant curvature scalar to be 




(2.5) 
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with Af the cosmo logical constant due to the f{R) gravity. Substituting this expression 
into (12.31) . one obtains the Ricci tensor 

The constant curvature black hole solution is given by 

dsl = - fl - ^ - ^Adt' + ■ + r^dnl, (2.7) 

r 3 

where > 0,< 0, =0 denote the Schwarzschild-de Sitter(dS), Schwarzschild-anti de Sit- 
ter (AdS), and Schwarzschild black holes, respectively. We call these "/(-R) black holes" 
because these were obtained from f{R) gravity. Even though these solutions are also ob- 
tained from the Einstein gravity with cosmological constant, their perturbation analysis is 
different from the Einstein gravity. The Schwarzschild-dS black hole has been extensively 
studied in f{R) gravities together with the cosmological implications of f{R) gravities [7]. 
However, as far as we know, there was a perturbation study on the Schwarzschild black 
hole in f{R,G) gravities |20j . 

Now we introduce the perturbation around the constant curvature black hole to study 
stability of the black hole 

9i.,y = 9ixv + h^,y. (2.8) 

Hereafter we denote the background quantities with the "overbar" . The linearized equation 
to (12.21) is given by (requiring f"{R) ^ 0) 



f{R)5R^^{h)-i^h,,+nR) 



SR{h) = 0,(2.9) 



2f"{Ry 

where the linearized Ricci tensor and curvature scalar take the forms 

SRf^uih) = i(^V^VAp + V^V,/i^p- V'V- V,.V,/i), (2.10) 
6R{h) = V^V" V - V^/i - A//i. (2.11) 

In order to find the black hole solution, we have to choose a specific form of f{R) as [18] 

f{R)=aiR + a2R^ + a3R^ + --- . (2.12) 
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We check that /(O) = at i? = -R = 0, which corresponds to either the Schwarzschild or 
Kerr black hole solution. In this work, for simplicity, we select a spherically symmetric 
Schwarzschild black hole solution with Aj = 0, 

dsl^^ = g^^dx^'dx" = -e^^'^de + e-^'^'^dr^ + r^{de^ + sin^ Od^p^) (2.13) 

with 

= 1 _ ^. (2.14) 

r 

Then, the linearized equation fl2.9p together with fl2.10p and (12.111) becomes 

+ [jF^] - V^V,) (v^V'^/i^^ - V^/i) = 0. (2.15) 

Taking the trace of (I2.15P with g'^", one has the fourth order equation for /i^^ 

- m^) (W'^V^K^ - V'/i) = (2.16) 
with the mass squared m'j defined by 

m} - (2.17) 

At this stage, we note that it is not easy to make a further progress on the perturbation 
analysis because there exist fourth order derivatives. We mention that for the Einstein 
gravity with f{R) = R, f'{R) = 1 and /"(O) = 0. In this case, one finds the equation 
for linearized curvature scalar: 6R{h) = 0, which means that 6R{h) is not a physically 
propagating mode. Actually, this equation leads to one constraint 

V^'V^h^p = V^h (2.18) 

which will also be recovered from the transverse gauge. Up to now, we did not fix any 
gauge. We would like to comment on the linearized equation when choosing the Lorentz 
gauge 

\/^h^'' = iv'^/i. (2.19) 
Under this gauge-fixing, the linearized equation fl2.15p takes the form [19] 



+ 2R^p,„hP'' + {g^^V^ - V.) V^h = (2.20) 



with the trace-reversed perturbation h^^ = h^^ — [21]. This equation was mainly 

used to mention that perturbed Kerr black holes obtained from f{R) gravity can probe 
deviations from the Einstein gravity [19]. Even though equation f l2.20p is simpler than 
f l2.15p . it is a non-trivial task to decouple odd and even perturbations around the Kerr 
black hole, arriving at two fourth order equations hopefully. Furthermore, we do not know 
how to solve the fourth order differential equation. Finally, we may choose the transverse 
gauge which works well for studying the graviton propagations on the the AdS4 spacetime 
background [221 ES] 

V^/i/^'^ = V'^/i, (2.21) 

which leads to fl2.18p when operating V on both sides. Using the relation f l2.18p . one 
immediately finds that the effect of f{R) gravity [2/"(0)-term in f l2.15p ] disappears because 
of 6R{h) = 0, leading to the Einstein gravity. Hence, the non-covariant gauge-fixing for the 
black hole perturbation should be different from the covariant gauge-fixing for the graviton 
propagations on the AdS, dS and Minkowski spacetimes [21]. It seems that the best way to 
resolve the difficulty confronting with the fourth order equation is to translate the fourth 
order equation into the second order equations by introducing auxiliary scalar fields. In 
other words, we must make a transformation from /(-R) gravity to the scalar-tensor theory 
(like Brans-Dicke theory) to analyze the stability of f{R) black hole. 



3 Perturbation of the scalar-tensor theory 

In this section, we will develop the perturbation analysis around the f{R) black holes (12. 7p 
in the different frame, the scalar-tensor theory. Introducing two auxiliary fields and A, 
one can rewrite the action (12. ip as [211 123 [2] 

Sst = i^,j d'x^g {0 {R-A) + f{A)} . (3.1) 

Varying for the fields and A lead to two equations 

R = A, <l> = f'iA). (3.2) 

Note that imposing (13.21) on the action (13.11) recovers the original action (12. ip . On the other 
hand, the equation of motion for the metric tensor can be obtained by 

(f)R^, - ^9^, + (^g^uV^ - V,) = 0. (3.3) 



In deriving the above equation, we used R = A whose rehabihty was discussed to explain 
the solar system test of /(-R) garvity [2]. In this work, we use this relation to analyze the 
stability of f{R) black hole only. Considering a constant curvature scalar R = R = A 
together with = f'{A) = const, Eq. fl3.3p becomes 



f{A)R^, - -g,JiA) = 0. 



Taking the trace of (13 ■4p leads to 



Rf'{A) - 2f{A) = 



(3.4) 



(3.5) 



which determines the positive, negative and zero curvature scalar by choosing a form of 



R 



4A4. 



(3.6) 



Substituting this expression into (13. 4p . one finds the Ricci tensor which determines the 
maximally symmetric Einstein spaces including Minkowski space 

1 /(^) . 



R,, 



(3.7) 



Now we are in a position to study the perturbation around the constant curvature black 
hole (EH). In addition to ([231), from (O, we have 



R + 6R{h) = A + 6A, 



f\A) + r\A)5A, 



(3.8) 



which leads to 

5R{h) ^ 5 A, 5(j) f"{A)5A. (3.9) 

Thus, instead of SR{h) and 50, we use SA as a perturbed field in addition to h^u. We 
expect that the same results can be derived when using /i^,^ and 50 in the Brans-Dicke 
theory because 50 ^ 6A. 

The linearized equation to (13. 3p takes the form 



f"{A)f{A)-f\A) 
2f'\A) 



5A 



m) J 
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(3.10) 



where the hnearized Ricci tensor 5R^y{h) is given by (12 .101) . It is important to note that 
taking the trace of fl3.10p with g^^'^ leads to the hnearized second order "scalaron" equation, 
instead of the hnearized fourth order curvature scalar equation (12.161) . as 



ml = 0, 



where the scalaron mass squared is given by 

r\A)-2f{A)nA) f'{A) 



3A., 



(3.11) 



(3.12) 



3/'(A)/"(A) 3/"(A) 

which was already known as (97) of [2] in dS spacetimes. This is the main why we have 
introduced the action (13. ip . instead of (12. ip . Plugging Eq. (l3.11|) into Eq. (l3.10p and rear- 
ranging the terms, we arrive at the linearized second order equation 



f'^jA) + fiA)f"{A) 
6/'2(A) 



5A. 



At this stage, we introduce a specific form of f{A) inspired from (I2.12p [T8 

f{A) = aiA + a2A^ + 03^^ ■ ■ ■ , 



(3.13) 



(3.14) 



where oi, 02, CI3, ■ ■ • are arbitrary constants. Their mass dimensions are [ai] = 0, [02] = —2 
because [/(^)] = 2, [A] = 2, and [0] = 2. In this work, we confine ourselves to the 
asymptotically flat spacetimes with = which accommodates the Schwarzschild black 
hole. In this case, we have to choose A = and thus. 



/(O) = 0, /'(O) = ai, no) = 2a2. 



Taking into account (I3.15p . Eq. (l3.13p reduces to 

1 



3m\ 



1 



where the scalaron mass squared in the asymptotically fiat spacetimes is given by 

f(0) 



(3.15) 



(3.16) 



(3.17) 



3/"(0)' 

Since the mass dimension of the linearized scalaron is two {[SA] = 2), it would be better to 
write the canonical linearized equations by introducing a dimensionless scalaron 6 A defined 
by 



6A 



5A 
3m\' 



(3.18) 



Finally, we arrive at two linearized equations 



V'-m\J5A = 0, 



m, 



g^JA = 0, 



(3.19) 
(3.20) 



which are our main result for carrying out the stability analysis of f{R) black hole. Impor- 
tantly, we observe that when replacing 



SA — >■ ip, m\ — )■ 



(3.2i; 



Eqs.f l3.19|) and f l3.20p are exactly the same equations of the massive Brans-Dicke theory for 
the stability analysis of the Schwarzschild black hole [26] . We mention that the stability 
analysis for the Schwarzschild black hole in the Brans-Dicke theory without potential (c = 0) 
has been established in [27j. Hence, we will use the result for the stability analysis for the 
Schwarzschild black hole in the massive Brans-Dicke theory (c 7^ 0). 



4 Stability analysis of f[R) black hole 
in the scalar-tensor theory 

The metric perturbations h^j,^ are classified depending on the transformation properties 
under parity, namely odd (axial) and even (polar). Using the Regge- Wheeler [28], and 
Zerilli gauge [21] , one obtains two distinct perturbations : odd and even perturbations. 
For odd parity, one has with two off-diagonal components and hi 

/ /io(r) \ 

hi{r) 



V ho{r) hi{r) ) 

while for even parity, the metric tensor takes the form with four components i^O) Hi, H2, 
and K as 

/ Ho{r)e''^''^ Hi{r) \ 

Hi{r) H2{r)e-'''^'"> 
PK{r) 

\ r'^sin^OKir) J 



-ikt 



sin 9 



dpi 



(4.1) 



(4.2) 
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where pi is Legendre polynomial with angular momentum / and e'^^''^ was given by fl2.14p . 
We note that fl4.ip and fl4.2p correspond to the non-covariant gauge- fixing [5U] . 

In order to explain the stability analysis of f{R) black hole briefiy, we mainly use the 
result for the stability analysis for the massive Brans-Dicke theory [26]. For the odd-parity 
perturbation, its linearized equation takes a simple form as 



5R^,{h) = 



(4.3) 



which leads to the Regge- Wheeler equation by introducing the tortoise coordinate r* 
r + 2mln[r /2m — 1] 



dr 



*2 



Q = 0, 



(4.4) 



6m 

.3 



(4.5) 



where the Regge- Wheeler potential is given by [28] 

VRw{r) = I 1 ) — 

This potential is always positive for whole range of — oo < r* < oo and a barrier- type 
localized around r* = 0, (see Fig. 3 of [27]) which implies that the odd-perturbation is 
stable [H]. For the even-perturbation, we have to use the linearized equation f l3.20p because 
the scalaron 6 A contributes to making an even mode M together with Hq, Hi, H2, and K 



definitely. After a long algebraic manipulation, we arrive at the Zerilli's equation 



dr 



*2 



+ 



Vz 



M = 0, 



(4.6) 



where the Zerilli potential is given by [29] 



Vzir) 



2m\ 
r J 



with 



2X^{X + l)r^ + QX^mr'^ + ISXm^r + 18m^ 
r3(Ar + 3m)^ 



A = -(/-!)(/ + 2). 



(4.7) 



(4.^ 



The Zerilli potential Vz is always positive for whole range of — 00 < r* < 00 and a barrier- 
type localized around r* = (see Fig. 3 of [27]), which implies that the even-perturbation 
is stable, even though the scalaron 6A is coupled to the even-parity perturbations. 
Finally, when considering 



(4.9) 
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the linearized scalaron equation fl3.19p leads to the Schrodinger-type equation 



where the scalaron potential is given by 



^ = 0, 



VA(r) 



2m\ 
r J 



(4.10) 



(4.11) 



where the second term is the usual scalar term with spin zero [in general, — 2m(s^ — 1) 
for s spin- weight of the perturbing field] , while the last term shows clearly the feature of a 
massive scalaron arisen from f{R) gravity. The potential Va{t) is always positive exterior 
the event horizon if the mass squared m\ is positive (non-tachyonic mass). For the shapes 
of scalaron potential, see Fig. 2 and Fig. 3 of by replacing m\ by — c. However, if m\ is 



negative (tachyonic mass), the potential becomes negative for large r and approaches — 

at infinity, indicating that the f{R) black hole is unstable against the scalaron-perturbation. 



5 Discussions 

We have investigated the stability of f{R) (Schwarzschild) black hole obtained from the 
f{R) gravity. Actually, it seems to be a formidable task to carry out the perturbation 
analysis around the black hole because the linearized Einstein equation is fourth order in 
f{R) gravity. 

We have proposed that the best way to resolve the difficulty confronting with the fourth 
order differential equation is to translate the fourth order equation into the second order 
equation by introducing auxiliary scalar fields. In this case, the linearized curvature scalar 
SR{h) becomes a massive scalaron, showing that all linearized equations are second order. 
We observed that the canonical linearized equations become the same equations for the 
massive Brans-Dicke theory when replacing the scalaron (its mass m\) by the Brans-Dicke 
scalar (its mass — c). Then, it is straightforward to make a decision on the stability of f{R) 
black hole. 

The stability on the metric perturbations remains unchanged, confirming that the odd 
(even) perturbations lead to the Schrodinger-type equation with the Regge- Wheeler (Zerilli) 
potential. This corresponds to the Einstein gravity, even though the even mode contains 
the scalaron in addition to Hq, Hi, H2, and K. The difference comes from the linearized 
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scalaron equation because the scalaron is a massive scalar which is physically propagating 
on the black hole background. It turns out that the f{R) is stable against the external 
perturbations if the scalaron does not have a tachyonic mass (/"(O) > 0). This is consis- 
tent with other perturbation analysis: the Dolgov-Kawasaki instability with f"{R) < in 
cosmological perturbations |3T], graviton and scalar propagations in the Minkowski [21], 
dS [521 ES] and AdS [23] spacetimes. 

In this work, even though the perturbation formalism is suitable for all constant curva- 
ture black holes (12 .Tp . including the Kerr black hole, we have analyzed the Schwarzschild 
black hole. We conjecture that the stability of black holes from Einstein gravity theory may 
hold for the f{R) black holes if one uses the scalar-tensor approach developed in Section 3. 
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1 Introduction 



Modified gravity theories, f{R) gravities [H El [3] have much attentions as one of strong 
candidates for explaining the current accelerating universe f{R) gravities can be con- 
sidered as Einstein gravity (massless graviton) with an additional scalar. For example, it 
was shown that the metric-/(i?) gravity is equivalent to the wbd = Brans-Dicke (BD) 
theory with the potential [5]. Although the equivalence principle test in the solar system 
imposes a strong constraint on f{R) gravities, they may not be automatically ruled out if 
the Chameleon mechanism is introduced to resolve it. It was shown that the equivalence 
principle test allows f{R) gravity models that are indistinguishable from the ACDM model 
in the background universe evolution [6]. However, this does not necessarily imply that 
there is no difference in the dynamics of perturbations [7]. 

In order for f{R) gravities to be acceptable, they must obey certain minimal require- 
ments for theoretical viability [H [S] . Three important requirements are included: (i) they 
possess the correct cosmological dynamics, (ii) they are free from instabilities (tachyon) 
and ghosts [8l|9l[6], (iii) they attain the correct Newtonian and post-Newtonian limits. 

On the other hand, the Schwarzschild-de Sitter black hole was obtained for a positively 
constant curvature scalar in [7] and other black hole solution was recently found for a 
non-constant curvature scalar [10]. A black hole solution was obtained from f{R) gravities 
by requiring the negative constant curvature scalar R = R [H]. If 1 -|- f'{R) > 0, this 
black hole is similar to the Schwarzschild-AdS (SAdS) black hole. In order to obtain the 
constant curvature black hole solution from f{R) gravity coupled to the matter, the trace 
of its stress-energy tensor T^,y should be zero. Hence, two known matter fields are the 
Maxwell [I2] and Yang-Mills fields [I3]. 

All black hole solution must pass the stability test. A black hole solution should be stable 
against the external perturbations because it stands as the physically realistic object [2]. 
There are two ways to achieve the stability of a black hole: one is the full stability by 
considering odd and even perturbations [15] and the other is the restricted stability by 
taking into account the spherically symmetric perturbations for simplicity [16]. The latter 
is not enough to guarantee the full stability and thus, it must be supported by the black 
hole thermodynamics (heat capacity) [1^. The basic idea is to decouple the linearized 
(perturbed) equations and then, manage to arrive at the second order Schrodinger-type 
equations for the physical field with the potential. If all potentials are positive for whole 
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range outside the event horizon, the black hole under the consideration is stable. Studies of 
stability of Kerr black hole are not as straightforward [IB] , because it is axially symmetric 
black hole and thus, the decoupling process seems to be complicated. However, this method 
is not suitable for f{R) black holes because f{R) gravity contains fourth order derivatives 
in the linearized equations [191 [20]. In this case, the requirement (ii) will play an important 
role in testing the stability of f{R) black holes. 

In this work, we investigate the stability of f{R) (Schwarzschild) black hole. We trans- 
form /(-R) gravity into the scalar-tensor theory to eliminate fourth order derivative terms by 
introducing two auxiliary scalars. Then, the linearized curvature scalar becomes a scalaron, 
indicating that all linearized equations are second order. Interestingly, they are exactly the 
same equations for the massive Brans-Dicke theory. Using the stability analysis of black 
hole in the massive Brans-Dicke theory, we show clearly that the f{R) black hole is stable 
against the external perturbations if the scalaron does not have a tachyonic mass. 

2 Perturbation of f(R) black holes 

Let us first consider f{R) gravity without any matter fields whose action is given by 



where ' denotes the differentiation with respect to its argument. It is well-known that 
Eq. (12.21) has a solution with constant curvature scalar R = R. In this case, Eq. (12. 2p can 
be written as 




(2.1) 



where k"^ 



SttG. The Einstein equation takes the form 




(2.2) 



RfiufiR) 



^g,J{R) = 0, 



(2.3) 



and thus, the trace of (12. 3p becomes 



RfiR) 



2f{R) = 0. 



(2.4) 



Note that the above equation determines the constant curvature scalar to be 




(2.5) 
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with Af the cosmo logical constant due to the f{R) gravity. Substituting this expression 
into (12.31) . one obtains the Ricci tensor 

The constant curvature black hole solution is given by 

dsl = - fl - ^ - ^Adt' + ■ + r^dnl, (2.7) 

r 3 

where > 0,< 0, =0 denote the Schwarzschild-de Sitter(dS), Schwarzschild-anti de Sit- 
ter (AdS), and Schwarzschild black holes, respectively. We call these "/(-R) black holes" 
because these were obtained from f{R) gravity. Even though these solutions are also ob- 
tained from the Einstein gravity with cosmological constant, their perturbation analysis is 
different from the Einstein gravity. The Schwarzschild-dS black hole has been extensively 
studied in f{R) gravities together with the cosmological implications of f{R) gravities [7]. 
However, as far as we know, there was a perturbation study on the Schwarzschild black 
hole in f{R,G) gravities |21j . 

Now we introduce the perturbation around the constant curvature black hole to study 
stability of the black hole 

9i.,y = 9ixv + h^,y. (2.8) 

Hereafter we denote the background quantities with the "overbar" . The linearized equation 
to (12.21) is given by (requiring f"{R) ^ 0) 



f{R)5R^^{h)-i^h,,+nR) 



SR{h) = 0,(2.9) 



2f"{Ry 

where the linearized Ricci tensor and curvature scalar take the forms 

SRf^uih) = i(^V^VAp + V^V,/i^p- V'V- V,.V,/i), (2.10) 
6R{h) = V^V" V - V^/i - A//i. (2.11) 

In order to find the black hole solution, we have to choose a specific form of f{R) as [19] 

f{R)=aiR + a2R^ + a3R^ + --- . (2.12) 
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We check that /(O) = at i? = -R = 0, which corresponds to either the Schwarzschild or 
Kerr black hole solution. In this work, for simplicity, we select a spherically symmetric 
Schwarzschild black hole solution with Aj = 0, 

dsl^^ = g^^dx^'dx" = -e^^'^de + e-^'^'^dr^ + r^{de^ + sin^ Od^p^) (2.13) 

with 

= 1 _ ^. (2.14) 

r 

Then, the linearized equation fl2.9p together with fl2.10p and (12.111) becomes 

+ [jF^] - V^V,) (v^V'^/i^^ - V^/i) = 0. (2.15) 

Taking the trace of (I2.15P with g'^", one has the fourth order equation for /i^^ 

- m^) (W'^V^K^ - V'/i) = (2.16) 
with the mass squared m'j defined by 

m} - (2.17) 

At this stage, we note that it is not easy to make a further progress on the perturbation 
analysis because there exist fourth order derivatives. We mention that for the Einstein 
gravity with f{R) = R, f'{R) = 1 and /"(O) = 0. In this case, one finds the equation 
for linearized curvature scalar: 6R{h) = 0, which means that 6R{h) is not a physically 
propagating mode. Actually, this equation leads to one constraint 

V^'V^h^p = V^h (2.18) 

which will also be recovered from the transverse gauge. Up to now, we did not fix any 
gauge. We would like to comment on the linearized equation when choosing the Lorentz 
gauge 

\/^h^'' = iv'^/i. (2.19) 
Under this gauge-fixing, the linearized equation fl2.15p takes the form [20] 



+ 2R^p,„hP'' + {g^^V^ - V.) V^h = (2.20) 



with the trace-reversed perturbation h^^ = h^^ — [22]. This equation was mainly 

used to mention that perturbed Kerr black holes obtained from f{R) gravity can probe 
deviations from the Einstein gravity pO]. Even though equation f l2.20p is simpler than 
f l2.15p . it is a non-trivial task to decouple odd and even perturbations around the Kerr 
black hole, arriving at two fourth order equations hopefully. Furthermore, we do not know 
how to solve the fourth order differential equation. Finally, we may choose the transverse 
gauge which works well for studying the graviton propagations on the the AdS4 spacetime 
background [231 El] 

V^/i/^'^ = V'^/i, (2.21) 

which leads to fl2.18p when operating V on both sides. Using the relation f l2.18p . one 
immediately finds that the effect of f{R) gravity [2/"(0)-term in f l2.15p ] disappears because 
of 6R{h) = 0, leading to the Einstein gravity. Hence, the non-covariant gauge-fixing for the 
black hole perturbation should be different from the covariant gauge-fixing for the graviton 
propagations on the AdS, dS and Minkowski spacetimes [22]. It seems that the best way to 
resolve the difficulty confronting with the fourth order equation is to translate the fourth 
order equation into the second order equations by introducing auxiliary scalar fields. In 
other words, we must make a transformation from /(-R) gravity to the scalar-tensor theory 
(like Brans-Dicke theory) to analyze the stability of f{R) black hole. 



3 Perturbation of the scalar-tensor theory 

In this section, we will develop the perturbation analysis around the f{R) black holes ((27 
in the different frame, the scalar-tensor theory. Introducing two auxiliary fields and A, 
one can rewrite the action (12. ip as [23 [221 [2] 

Sst = i^,j d'x^g {0 {R-A) + f{A)} . (3.1) 

Varying for the fields and A lead to two equations 

R = A, <l> = f'iA). (3.2) 

Note that imposing (13.21) on the action (13.11) recovers the original action (12. ip . On the other 
hand, the equation of motion for the metric tensor can be obtained by 

(f)R^, - ^9^, + (^g^uV^ - V,) = 0. (3.3) 



In deriving the above equation, we used R = A whose rehabihty was discussed to explain 
the solar system test of /(-R) garvity [2]. In this work, we use this relation to analyze the 
stability of f{R) black hole only. Considering a constant curvature scalar R = R = A 
together with = f'{A) = const, Eq. fl5.3p becomes 



f{A)R^, - -g,JiA) = 0. 



Taking the trace of (13 ■4p leads to 



Rf'{A) - 2f{A) = 



(3.4) 



(3.5) 



which determines the positive, negative and zero curvature scalar by choosing a form of 



R 



4A4. 



(3.6) 



Substituting this expression into (13. 4p . one finds the Ricci tensor which determines the 
maximally symmetric Einstein spaces including Minkowski space 

1 /(^) . 



R,, 



(3.7) 



Now we are in a position to study the perturbation around the constant curvature black 
hole (EH). In addition to ([231), from (O, we have 



R + 6R{h) = A + 6A, 



f\A) + r\A)5A, 



(3.8) 



which leads to 

5R{h) ^ 5 A, 5(j) f"{A)5A. (3.9) 

Thus, instead of SR{h) and 50, we use SA as a perturbed field in addition to h^u. We 
expect that the same results can be derived when using /i^,^ and 50 in the Brans-Dicke 
theory because 50 ^ 6A. 

The linearized equation to (15. 3p takes the form 



f"{A)f{A)-f\A) 
2f'\A) 



5A 



m) J 
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(3.10) 



where the hnearized Ricci tensor 6R^y{h) is given by (12 .101) . It is important to note that 
taking the trace of fl3.10p with g^^'^ leads to the hnearized second order "scalaron" equation, 
instead of the hnearized fourth order curvature scalar equation (12.161) . as 

(3.11) 



V -mljM = 0, 
where the scalaron mass squared m\ is given by 

r\A)-2f{A)nA) f'iA) 



m. 



(3.12) 



3f(A)/"(A) Sf'iA) 

which was already known as (97) of [2] in dS spacetimes. This is the main reason why we 
have introduced the action (13. ip . instead of (12. ip . Plugging Eq. (l3.1ip into Eq. (l3.10p and 
rearranging the terms, we arrive at the linearized second order equation 



5Rf,^{h) - AaH 



r\A) + f{A)f"{A) 
6/'2(A) 



5A. 



(3.13) 



In this work, we confine ourselves to the asymptotically flat spacetimes with = which 
accommodates the Schwarzschild black hole. Also, we do not choose an explicit form of 
f{A) since such a restriction seems unnecessary to study the stability of f{R) black holes. 
In this case, we have to choose A = and thus, 



/(o) = o, /'(o)^o, r(o)^o. 



Taking into account (I3.14p . Eq. (l3.13p reduces to 

1 



3m\ 



Vf,VJA + -gf,JA, 



(3.14) 



(3.15) 



where the scalaron mass squared in the asymptotically flat spacetimes is given by 



(3.16) 



3/"(0)' 

Since the mass dimension of the linearized scalaron is two {[6 A] = 2), it would be better to 
write the canonical linearized equations by introducing a dimensionless scalaron 6 A defined 
by 



SA 



SA 
3m\' 



(3.17) 
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Finally, we arrive at two linearized equations 




(3.18) 




gf,JA = 0, 



(3.19) 



which are our main result for carrying out the stability analysis of f{R) black hole. Impor- 
tantly, we observe that when replacing 



Eqs.f l3.18|l and f l3.19p are exactly the same equations of the massive Brans-Dicke theory for 
the stability analysis of the Schwarzschild black hole [27]. We mention that the stability 
analysis for the Schwarzschild black hole in the Brans-Dicke theory without potential (c = 0) 
has been established in . Recently, the scalar field perturbations of Schwarzschild black 
hole was carried out in tensor- vector-scalar theory [29]. Hence, we will use the result for 
the stability analysis for the Schwarzschild black hole in the massive Brans-Dicke theory 



Finally, we would like to mention that even though the Brans-Dicke scalar was intro- 
duced instead of scalaron 6 A, there is no change in the linearized equations. (See Appendix). 
For this purpose, we rewrite fl3.13p as 



The only difference seems to be the mass squared m^, in compared to m\ in (13.121) . However, 
they become the same mass squared when making replacements (I5.14p . 

4 Stability analysis of f{R) black hole 
in the scalar-tensor theory 

The metric perturbations h^^ are classified depending on the transformation properties 
under parity, namely odd (axial) and even (polar). Using the Regge- Wheeler [5U], and 
Zerilli gauge [21] , one obtains two distinct perturbations : odd and even perturbations. 



6A — )■ Lf, m\ —f — c. 



(3.20) 



(c^O). 



(3.21) 
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(4.1) 



For odd parity, one has with two off-diagonal components /iq and hi 

/ /io(r) \ 

hi{r) 



V h^{r) h{r) / 

while for even parity, the metric tensor takes the form with four components Hq, Hi, H2, 
and K as 



^^*sin^^ 
dO 



Hi{r) 
























(4.2) 



sin^ ^Js:(r) / 

where pi is Legendre polynomial with angular momentum / and e^^"^^ was given by fl2.14p . 



We note that (14.1 p and (14. 2 p correspond to the non-covariant gauge-fixing 

In order to explain the stability analysis of f{R) black hole briefly, we mainly use the 
result for the stability analysis for the massive Brans-Dicke theory [27] . For the odd-parity 
perturbation, its linearized equation takes a simple form as 



8R^,{h) = 



(4.3) 



which leads to the Regge- Wheeler equation by introducing the tortoise coordinate r* 
r -|- 2m ln[r/2m — 1] 



dr 



*2 



+ 



RW 



Q = 0, 



(4.4) 



6m 

.3 



(4.5) 



where the Regge- Wheeler potential is given by [30] 

VRwir) = I 1 ) — 

This potential is always positive for whole range of —00 < r* < 00 and a barrier-type 
localized around r* = 0, (see Fig. 3 of [28]) which implies that the odd-perturbation is 
stable [in]. For the even-perturbation, we have to use the linearized equation (I3.19P because 
the scalaron 6 A contributes to making an even mode M together with i^O) Hi, H2, and K 
definitely. After a long algebraic manipulation, we arrive at the Zerilli's equation 

d'^M 



dr 



*2 



k'-Vz 
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M = 0, 



(4.6) 



where the Zerilh potential is given by [31] 

2A2(A + l)r^ + + ISAmV + ISm^ 



Vz r 



2m \ 
r / 



r^(Ar + 3m)" 



with 



A = -(/-l)(/ + 2). 



(4.7) 



(4. 



The Zerilh potential Vz is always positive for whole range of — cxd < r* < cxd and a barrier- 
type localized around r* = (see Fig. 3 of [2B]), which implies that the even-perturbation 
is stable, even though the scalaron 6A is coupled to the even-parity perturbations. 
Finally, when considering 



5AocJ:^Yim{e,ip)e-'''\ 
r 

the linearized scalaron equation fl3.18p leads to the Schrodinger-type equation 



+ 



^ = 0, 



where the scalaron potential is given by 



VAir) =1 



2m\ 



/(/ + i) 



2m 



m, 



(4.9) 



(4.10) 



(4.11) 



where the second term is the usual scalar term with spin zero [in general, — 2m(s^ — 1) 
for s spin- weight of the perturbing field] , while the last term shows clearly the feature of a 
massive scalaron arisen from f{R) gravity. The potential VA{r) is always positive exterior 
the event horizon if the mass squared m\ is positive (non-tachyonic mass). For the shapes 
of scalaron potential, see Fig. 2 and Fig. 3 of [27] by replacing m\ by — c. However, if m^ is 
negative (tachyonic mass), the potential becomes negative for large r and approaches — |m^| 
at infinity, indicating that the f{R) black hole is unstable against the scalaron-perturbation. 



5 Discussions 

We have investigated the stability of f{R) (Schwarzschild) black hole obtained from the 
f{R) gravity. Actually, it seems to be a formidable task to carry out the perturbation 
analysis around the black hole because the linearized Einstein equation is fourth order in 
f{R) gravity. 
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We have proposed that the best way to resolve the difficuhy confronting with the fourth 
order differential equation is to translate the fourth order equation into the second order 
equation by introducing auxiliary scalar fields. In this case, the linearized curvature scalar 
6R{h) becomes a massive scalaron, showing that all linearized equations are second order. 
We observed that the canonical linearized equations become the same equations for the 
massive Brans-Dicke theory when replacing the scalaron (its mass m^) by the Brans-Dicke 
scalar (its mass — c). Then, it is straightforward to make a decision on the stability of f{R) 
black hole. 

The stability on the metric perturbations remains unchanged, confirming that the odd 
(even) perturbations lead to the Schrodinger-type equation with the Regge- Wheeler (Zerilli) 
potential. This corresponds to the Einstein gravity, even though the even mode contains 
the scalaron in addition to Hq, Hi, H2, and K. The difference comes from the linearized 
scalaron equation because the scalaron is a massive scalar which is physically propagating 
on the black hole background. It turns out that the f{R) is stable against the external 
perturbations if the scalaron does not have a tachyonic mass (/"(O) > 0). This is consis- 
tent with other perturbation analysis: the Dolgov-Kawasaki instability with f"{R) < in 
cosmological perturbations [33], graviton and scalar propagations in the Minkowski [25] . 
dS [Ml [33 and AdS [21] spacetimes. 

We would like to compare our stability analysis with ref. [11], where the very restricted 
stability was performed by taking into account the spherically symmetric static pertur- 
bations for simplicity. Therefore, their results did not show the full stability analysis, in 
contrast with our results. 

In this work, even though the perturbation formalism is suitable for all constant curva- 
ture black holes (12. 7p . including the Kerr black hole, we have analyzed the Schwarzschild 
black hole. We conjecture that the stability of black holes from Einstein gravity theory may 
hold for the f{R) black holes if one uses the scalar-tensor approach developed in Section 3. 
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Appendix: Brans-Dicke approach 

We may rewrite ( 13. ip as the Brans-Dicke theory 

= ^ y d^x^g {<PR - Vm (5.1) 

with the potential 

y(0) = 0A(0)-/(A(0)). (5.2) 
Varying for the fields g^^^ lead to the following equations: 



(5.3) 



R = V'{(j)) (5.4) 

where ' denotes differentiation with respect to 0. Taking the trace of 05. 3p in order to 
replace R in (15. 4p leads to the scalar equation 

3V'^(f) + 2V{(f))-(f)V'{(j)) =0. (5.5) 

For the constant curvature scalar case only, we have 

= 0, R = v'{^) = '^^^ = 4A^, R^, = k^g^, (5.6) 

where 

A<A = -| (5.7) 
with £ the AdS4 curvature radius. From (15. 4p and (15.61) . the potential may take the form 

V^ = Co + Ci(0 - 0) + C2(0 - 0)2 + ■ ■ ■ , (5.8) 

where 

Co = |-0 < 0, ci < 0, C2 > 0, > 0. (5.9) 

The linearized equations around the constant curvature scalar background (15. 6p can be 
obtained as 

(v2-mj)50 = O, (5.10) 

2 

5R^,{h) - A<^V - V^V,50 - + K\9^.u^ = 0, (5.11) 
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where 50 = 50/0 and is given by 

ml = \ (0V"(0) - V'(0)) , (5.12) 

where 

0-l/"(0) = 2c20 > 0. (5.13) 

We emphasize that the mass squared fl5.12p is exactly the same as one derived in the 
hterature [251 ESj with zero Brans-Dicke parameter (cu = 0). Also, it is easily shown that 
when making the replacements 

50 ^ 5i, ^ /'(A), \/(0) ^ /(A), \/'(0) ^ A, \/"(0) ^ (5.14) 

f l3.18p and f l3.2ip lead to fl5.10p and (15.110 . respectively. Hence, it is enough to solve (I3.18P 
and (I3.19P for the stabihty analysis of /(i?)-Schwarzschild black holes with A^i = 0. 
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